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Abstract 


In this paper, a two-dimensional time-fractional telegraph equation is 
considered with derivative in the sense of Caputo and 1 < 8 < 2. The aim 
of this work is to extend the Crank—Nicolson method for this time-fractional 
telegraph equation. The stability and convergence of the numerical method 
are investigated. Also, the accuracy and efficiency of the proposed method 
are demonstrated by numerical experiments. 
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1 Introduction 


Fractional calculus can be used to model many complex problems. It has 
been used in many fields of science, engineering, and finance [1, 4, 18, 25, 26]; 
this fact is the main source of inspiration for most of the recent studies 
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conducted on fractional calculus. The classical telegraph equation is used 
in random walk theory [2]. The time-fractional telegraph equation (TFTE) 
models the neutron transport process in the core of a nuclear reactor [27, 28]. 


In recent decades, the fractional telegraph equation has been solved by 
many researchers. Here, we briefly describe some of the studies that have 
been conducted in this field of research. Chen, Liu, and Anh [5] proposed 
the analytic solution of the TFTE using the separating variables method. 
The homotopy analysis method was used for the TFTE by Das et al. [6]. 
Yildirim [31] applied the homotopy perturbation method to solve space- and 
time-fractional telegraph equations. Momani [17] used Adomian decompo- 
sition methods to obtain the analytic and approximate solutions of space- 
and time-fractional telegraph equations. Biazar, Ebrahimi, and Ayati [3] 
proposed the variational iteration method to solve the fractional telegraph 
equation. Jiang and Lin [11] presented the exact solution of the TFTE using 
the reproducing kernel theorem. Nikan, Avazzadeh, and Machado [19] used a 
mesh-free spectral approach based on LRBF-FD to solve the TFTE with the 
fractional derivative described in the sense of Caputo. A radial basis function 
collocation method was used for solving the nonlinear TFTE by Sepehrian 
and Shamohammadi [22]. Hosseini et al. [9, 10] considered the meshless local 
radial point interpolation method, and Mohebbi, Abbaszadeh, and Dehghan 
[16] used the radial basis function technique for the TFTE. Shivanian applied 
spectral meshless radial point interpolation methods in [23], and the meshless 
local Petrov-Galerkin scheme was used in [24] to approximate the TFTE. 


Many researchers have studied the fractional telegraph equation by using 
the finite difference method. Liang, Yao, and Wang [15] considered the TFTE 
by using a fast, high-order difference scheme. The finite difference method 
was used to solve the linear TFTE by Li and Cao [14]. Wang and Mei [29] 
considered the TFTE using a Legendre spectral Galerkin method in space and 
the generalized finite difference scheme in time. For a time-space-fractional 
telegraph equation, Zhao and Li [32] used a finite difference method in time 
and a Galerkin finite element method in space. A numerical method for the 
TFTE was proposed by Wei, Liu, and Sun [30], in which they discretized this 
equation with a new finite difference scheme in time and a local discontinuous 
Galerkin (LDG) method in space. 


In this work, we find an approximate solution to the following TFTE [13]: 


OF u(z,y,t) | 0°*ula,y,t) 
ate ote-1 


+ u(z, y,t) = Au(z, y,t) + f(x,y, t) 
(x,y) €ACR’,0<t<T, (1) 


with initial and boundary conditions 


u(x, y,0) = y(a,y), (z,y) €N=NUAQ, (2) 
Puls U-9) _ (x,y), (x,y) €%=NUAQ, (3) 
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u(x, y,t) = h(a,y,t), (x,y) € OQ,t > 0, (4) 


where 1 < 8 < 2,A is the Laplace operator, OQ is the boundary of Q, 
f(x,y), e(a, y), v(a,y), and h(a,y,t) are continuous functions, u(x,y,t) € 
C?(Q x [0,7]) is an unknown function, and the fractional derivatives are 
defined in the sense of Caputo, as follows: 


OP-lu(a, y,t 1 * Ou(a, y, = 

wey =a | ey Mie s)!-Pds, 1<8<2, (5) 
Pula, y,t) it : u(x, y, s) 1-8 

7B = Tas) | 72 (t—s)* "ds, 1<B<2. (6) 


The Crank-—Nicolson difference scheme can be used easily for space- 
fractional equations, but some manipulations are needed for time-fractional 
equations [12]. In [19, 13], a semi-discrete scheme based on the Crank-— 
Nicolson method was used to discretize the time-fractional equation. In this 
work, the discretization of time-fractional derivatives is similar to [12]. The 
general idea for proving stability and convergence is taken from [19], but our 
approach differs from that in the details. 

The remainder of this paper is organized as follows. In Section 2, the 
discretization of (1) is described. The stability and the convergence of the 
proposed method are proved in Sections 3 and 4, respectively. Section 5 is 
devoted to the numerical tests. Finally, the conclusion is given in Section 6. 


2 Discretization of the problem 


In this section, we explain the discretization of (1) by using the Crank— 
Nicolson difference scheme, such that the proposed difference schemes are 
uniquely solvable. 

Consider Az and Ay as the grid sizes in space for the finite difference 
scheme, where {(xi, yi), 0; =iAzv, yj =jAy,0<i<1,0<j< J;I,J €R} 
covers 2. Also, N is a positive integer, and the grid size in time for the finite 

T 


difference scheme is At = 7. Assume that u/',; is the value of u (xj, 4j,tn) - 


The following lemma provides suitable tools for the discretization of (1). 


Lemma 1. If g(t) € C?[0,T] and 1 < 6 < 2, then 
(a) 


[ders - 98a 
1-6 
-oB C b+ 5)? -(n—b— 5) | [a (te) — 9 (tea) 


+ O(At)?-8, k=1,2,...,N—1. 
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ta 
[oe ltn4 - 9) Pas 
t 


(Apr? 


= pee ly lin) ~ 9 tna)] + OC Aiye-?, on EN. 


Proof. The Taylor expansion allows us to write 


g'(s) = g (tk) 7—g (te-1) 1 s) ” (m) 


At Day lt g — (te-1 — 8)" 9" (n2)], 


m € (8, te) ,N2 © (te-1, 8). 


It is easy to show that 


te 
(t, — 8)" (tn—2 — s)'-#ds = un (At)*-*, w, ER, 


tk-1 


tk 
/ (th-1 — 8) (t,4 - s)'-8ds = w(At)*-8, w2 ER. 


te—-1 


Thus, 
tk 
[ltr 4 9) Pas 


tk-1 
= tr 
_g (te) ae (t,4 — 8)! Pas 


tes)" bees 1 — s)'""ds 


Sof (th-1 — 8) nme —s)'-8ds 
th iP, 


ca g( (At)? = ies 1\o 
a og [tb phe 


w(At)?-F, weER. 


This completes the proof of part (a). Part (b) can be proved in the same 
way. 


By defining b, = (s as (s a SH 1 Desay WB <2; 
it is easy to show that 
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n-1 

_ cae Tt MY ohne, (ee 
Suh Wh) (mb SPP = (B= BI?) 4 ae (uly 085) 
k=1 


n—-2 

i oe 

— b-nt, + S- (by—K—-1 = bn—x) ur, + (sa-8 = bi urs 1 
k=1 


By using part (b) of Lemma 3, the discretization of (5) at the grid point 
(x;,y;) and the time step (1 — 4) is as follows: 


oP "ula, yyt) 1-4 1 ‘a Ou(xi, Yj, 8) 
B) Seo 


ae ‘ 1-8 
ae-1 td TQ 7 ies a 


(Ag? 1, bg 8) 


~ T(3— 8B) * 52-8 [uj uss + O(At)?-*. 


By using Lemma 3 and relation (7), the discretization of (5) at the grid point 


(xi, yj) and the time step (n — 3) is as follows: 


OP eg ised 1. eee Ou(xi, yz, 8) Fi 
19) oa = vz) ’ t Ba 
ate-1 li r(2 = B) », in Ae ( n—-4 s) 8 


1 tok Ou (Xi, 7,8) 1B 
7 T(2— ) i Os (20 8) ds 
_ (At)8 ie. eee ; 
“TB—A) —bn—1u) 5 — Da lrne — byw ub, 


1 n-1 1 n 
—(sa-8 a bi urs + atts} 


POMP, wed <o<t=1,1 35 271, 


(9) 
In addition, similar to (8) and (9), and by using the relation 
du|* ur; —uf5* At 0?u 
= ; ; 497) ) k 2 1, th it ’ 
a |, At tp pe Miryam) m € (th—1, tk) 
(10) 
we obtain 
1-i _ 0 
du(x,y,t)| 7 _ (At)? Se (ie ae Ou + O(At)2-%, 
oe |,, (3-8) 28 At Ot |, 
(11) 
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Pula, Y; t) me 
are td 


At 1-8 r) 0 n—2 uk , uke 
= : n—-1 . S- (Diy ea bn—k) ! ! 
T(3 — £) at |; = At 
n-1 n—-2 n n-1 
( 1 by) aid mar 1 Uag 7 Yay 
22-8 At 22-8 At 
LOA. weve rai oe f=1, 
(12) 
Having the Taylor expansion in mind, we can write 
ae + ur 2 
u(@isYirty—4) = 2G + Of), PS; 
n>11l<i<I-1ll<j<J-l, 
Nats! tue, 
Au(ti,¥jstyy) = ——? + O(a) 
ol tt — Qui, oe Tae Tarai _ up + ur 
2 (Aa)? (Ay)? 
Meg TAU Tg | bgt 7 Mig TUE 
(Ax)? (Ay)? 
+ O(Az)? + O(Ay)? + O(At)?, 
n>11l<i<I-1ll<j<J-1. 

(14) 


Using the finite difference schemes (11), (8), (13), and (14), the discretization 
of (1) at the grid point (a;,y;) and the time step (1 — $) is as follows: 


(At)!-8 1 


(At)!-8 i 
T(s— A)” BA 


T'(3 — 8) ‘ 92-6 (U3, u} 5) 


0 
Wig 7 Uig _ Ow 


At ot 


tJ 


0 0) 0 0 0 0 
Le, oy 1 J Cig 7 246g FUR | Vigta 724i + Yi j—1 
5 (ups + Ui5) = 


2 (Ax)? (Ay)? 


1 1 1 1 1 1 
vita Tig Ting | Cag 7 Mag Hit 


| (Ax)? | (Ay)? 


+ fiz? + O(Ac)? + O(Ay)? + O(At)?-*, 
a oe ee eae a Bo 
(15) 
Using the finite difference schemes (12), (9), (13), and (14), the discretization 
of (1) at the grid point (x;,y;) and the time step (n — 4) can be written as 
follows: 
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= 0 n—-2 k k-1 
T(3 a B) ot UJ k=1 At 
n—-1 n—2 n n-1 
( 1 bi) Ui j = Us J, Ui ZZ Ui 
22-8 At 22-8 At 
n—2 
(At)i-8 ; k 
tA 6 ith Oe 
— ud ae 
I(s— 6) hol 
n—-1 n 
1 n-1 1 n | Min i Med 
— (s-8 — bi)uy, tasty} D) 
_ lL fut Quip ty | aie a a 
D (Az)? (Ay)? 
4 Mag — Ug Tg, Ue = OU Ts 
(Aa)? (Ay)? 


+f"? + O(Ac)? + O(Ay)? + O(At)?*, 


metic tai 7S 7 =1, 
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(16) 


Finally, rearranging (15) and (16) and neglecting the truncation errors, it 
is obvious that the coefficient matrix of the unknowns is strictly diagonally 
dominant and so, it is nonsingular [8]. Therefore, by neglecting the truncation 
errors in (15) and (16), the unknowns [u?,| (1 <i<I-1,1<j<J-1) 


ig 


can be obtained for n = 1 and n > 2, respectively. Hence, the proposed 


Crank—Nicolson scheme is uniquely solvable. 


3 Stability 


In this section, we study the stability of the proposed Crank—Nicolson scheme 
for (1) with initial and boundary conditions (2)—(4). To do so, we introduce 
the following spaces and recall some theorems and lemmas, which will be 


used hereafter. 


H'(Q) = {ve 1°(Q): Dve L°(Q)}, 
HO) ={0 6. O)? Dolpg = 0} 
H?(Q) = {v € L7(Q) : D*ve L7(Q), jal < 2}. 


Theorem 1 (The Cauchy—Schwarz inequality). [21 


If wand v are members of an inner product space 2 with inner product (-,-), 


then 


S |lullz2|lullze. 


th 


(u,v) | = [wae 
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Theorem 2. (Green’s theorem) [21] 
If Q is a boundary domain in R”, then 


U 


[ vureax = | v5eds— f Auvdr, foru€ H?(Q),v € H'(Q). 
Q an O Q 


n 


Theorem 3 (The Poincare-Friedrich inequality). [21] 
Let Q“ be a boundary domain in R”. Then, there exists a constant c, > 0 
such that 

l|e||7.2 < cp||Vulle, for all u € Hj(Q). 


Theorem 4 (The discrete Gronwall theorem). [20] 
Assume that k, is a nonnegative sequence and that the sequence ¢,, satisfies 
the following relations: 


go < Jo; 


WV 
— 


n-1 n—1 
on Sg0+ > pet >> keds, nr 
s=0 s=0 
If go > 0 and py, > 0 (for n > 0), then 
n—1 n—1 
mrs (wtdr)an(Ss), me 
s=0 s=0 


We state some useful relations in Lemmas 2 and 3. These are easy to 
prove. 


ull |r|] <  |ull? 1 |lv||?, for all u,v € Q, for all y € 


Lemma 2. It holds that ye 


R. 


Lemma 3. If b, = ee aes (eee a (s =1,2,...,1< 8 < 2), then 
bn < On-1 <0 < bn < by <1. 


Neglecting the truncation errors, equations (15) and (16) can be written 


as 
0 
tJ 


sk 
(Aut; +Au?;)+f5?, 1<i<I-Ll<j<J-l, 


T(3—B) ~ 2-8 |~ At Ot 


(At)? 1 1 0 1 
= T@—B) © P (us.5 — Mig) + 


(At)!-6 1 E _ Us Ou 


1 
2 


and 
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(At)i-4 duly hg — eG 
Pe ge en es 
n—-1 n-2 n n—-1 

ee 

22-B At 92-8 At 

(At)!-8 : n—-2 ‘ 

ae oa ee S- (bn—K—1 — bn—k) Ug 3 

13 ~~ B) k=1 


1 n-1 1 n tJ tJ 
—(3=8 _ by Juz; + ats} + =r = 


(18) 
respectively. Let ub; (1 <i<I—-1,1<j < J—1,n =1,2,...) be the 
approximate solution of (17) and (18) with respect to the round-off error, 
and let up; (l1<i< I-11 <J—1,n=1,2,...) be the exact solution 
of (17) and (18). Define 
(O<i<I, 0<j<J, and n=0,1,...). 


a ~n 
Fe ae 


ae 


By considering e” instead of e7,, 


equations: 


Cues {{o* se + (e! o “(et te) = S(Ae! + Ac), 


we obtain the following round-off error 


T(3— 8) ~ 22-8 || At 
(19) 
(At)? a eae} 
T(3— 8B) —bn_10e — », (bn—h—-1 bn—k) At 
1 etl — er—2 1 et er—1 
(Sa-p ~ bs) Rp + = A \ 
(20) 


At)!-8 _ 1 
t ( ) bn—1e° oo (On—k—-1 = bn—k) e* — (sa-8 = by)e"—+ 


k=1 


1 x etl +e” Ae”—1 + Ae” 
8° \ 5 = 5 ; n> 2, 
0 7 )0 
where de? = ale sre oa ‘ye Now, we are ready to present the following 


theorem. 


Theorem 5. If e* € Hd(Q), then the solutions of the finite difference ap- 
proaches (17) and (18) are unconditionally stable. 


Proof. Let a= Oo: If we multiply (19) by (e’— e°), then we obtain 
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a a 1 
3-B(AD) (e —e°,e1 —e°) + a (eh — eel) + : (et + 69, el — e°) 
= ; (Act By Ae®,e" a gy) = 2 (de°, e* — a) : 


(21) 
Applying Theorem 2 (Green’s theorem) to (Ae! + Ae°,e! — e®) in the left 
side of (51) and applying Theorem 1 (the Cauchy—Schwarz inequality) and 
Lemma 2 to the right side of (51), we can write 


a a 1 
aap le - 2? + salle — 21 + sdle P= le) 
aca oy < of [de | let —e||" 
+ 5(llvell -Ivell) s | 2 1 2 
Therefore, 
2 2 2 a 2 
Iver <llel + Vell + sz [lee 


and by applying Theorem 3 (the Poincaré—Friedrich inequality), we find a 
constant c, > 0 such that 


QCp 
Oost 


|| Ver]? ~ (cp +1) ||Ve|)? + ||Voe ||”. (22) 


If we multiply (20) by (e” — e"~*), then we find 


a 


(e” ane meg _ er") 4 rae 


a n on n— 
22-B (At) oe 2 
+ ; (e? +e"-1 ce” — e” 1) 5 (Ae" + Ae"! e” — e” 1) 

=abn_} (5e°, e” — a 


n—-2 ek — ek-1 : 
Tra > (On—h—1 bn—k) ( ’ ef a ) 
At 
k=1 
1 etl _ er—2 


a( 


pep ~ 2) 


n-2 
+a S- (bn—k—-1 — bn—k) eben = en) 
k=1 


+ a( by) (e"~*,e" — e™*), 


1 
22-8 
(23) 
Applying Theorem 2 (Green’s theorem) to (Ae” + Ae™1,e™ — et} in the 
left side of (23) and applying Theorem 1 (the Cauchy—Schwarz inequality) 
and Lemma 2 to the right side of (23), we obtain 
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n n— a n a n n— 
le” =e MP (SP lien? — 5 (ere) 


oe 
2-F(At) 


i n n— 
+ 5(le"I? — |e") + 


1 
5 (lVer|? — ||ver™]]") 


< aby a(2 [seo]? + ee aan 
pes n—-1 2 272 
n-2 2 Lk k-1 
er Ie, 1 Ae pe ii2 
ee — bn—r) ( 5 | me \|? + a2 lle" — e”* ||") 
7 ters 24 n n—1))2 
of IP + lle” - 21) 
2 
7 2 e” —enrl 
+ abn_1( 5 lle || | rp | ) 
= ae k\|2 i n n—1||2 
aS acter dace) MAN ae PY 
fa(E er + aller — eH) yER 
2 27? ; ; 
(24) 
Furthermore, from Lemma 3, we deduce that 
ab as a . 2a 
n—-1 
gE (he hint hs GO eee, R. 25 
42 2 24 ( k-1 k) 2 7 YE ( ) 
Having (25) in mind, equation (23) gives 
sane ee aie | ea spas lle IP + 5lle I? + sllve I? 
¥? 2 2 ay? — 2 
cass {BI + PY +P Se aka bs) Hel 
4 OT nat —to-t) fet et? + Meg fer — oa 
2(At)? _ n—-k-1 n—k T 2(At)? 
ay? u n—1||2 2a n n—1}}2 a n on- 
tip t gh le sa le ie haa eve” 2) 
+5llve?, ver. 
(26) 
By using Theorem 1 (the Cauchy—Schwarz inequality) and Lemma 2, we 
obtain ; 
a on do a "||" Ce 
aie < sts! mr en \. 07) 
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Consider the following relations: 


Q2 n-2 
ay 442 ay en en 2 
Daa pe 2 (On-t-1 — Pn-x) et — ef ‘I + scape ile” * - eI 
Q n-2 
< a 5 Yo On—a-1 = bn) (Ile*ll? + fle T1)7) (28) 


1 
+ Tap (lle" IP + lle" 17) - 


If we use (27)-(28) and assume that y? = 2°-9(At), then relation (26) 
allows us to write 


1 a efe nil Wea nie 
5X sep llen? + 5 len? + 5 Ver 


¥ 2 2 ary? 2 
< abu {Hel + felt} + FS Gnade 


reid ae 2 14/2 ary? e- 
Tape Da On-ka — os) + le" |?) + err |" 
k=1 
bi ay? a oat 
("+ aap t ae t all alee ve tees 
(29) 


By using Theorem 3 (the Poincaré—Friedrich inequality), we find a constant 

Cp > 0 such that relation (29) implies 

abn—1y" wn 177 
2 


1 pyenpe <Q? a Ive? + 


2 n—2 
+ Joey 32 (bn—k—1 — bat) [Vel 
k=1 


“ 
apt n—k—-1 — bn—z)(||Verl] + || Vet )°) (30) 
2 


Gp || Ve" 


ay n—-2 ile ay? 
(At )2 2 op ae | + ( 9) 5s aia 
a 


+ 5x98 + 5)¢ || Ver ilies 5 lve 


+ 


n> 2. 


We may assume without loss of generality that there exist constants 01,02 > 0 
such that relations (22) and (30) can be written as 
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Ver] <A Ver’ + G2 || Voe"||”, 
n-1 
Ver? < (04 [fVe%]? + 6 [| 75?) +> (ce Vel) (3 
k=1 


n>2, 01,02>0, co >0 for k=1,...,n—-1. 


By Theorem 4 (the discrete Gronwall theorem), equation (31) yields 
5 s n—-1 
Ver ||? < (61 I|Ve>|| + 0 || Vse"|| ) exp (s +) , n>1,61,62 > 0, 
k=1 


and according to Theorem 3 (the Poincare-Friedrich inequality), there exists 
a constant ¢, > 0 such that 


lle"||? < i (61 Ive]? + 62 75°”) exp (© +) » n>1,01,02 >0. 
_ (32) 
By using Lemma 3, it is easy to show that 
oe >» Say? a 
dS (207 t Ape + ae 4 1) eth, (33) 


Set 0 = (207 ne + ae 4 1) Cp +1. Then, it follows from relations (32) 
and (33) that 


leh < fe (Au le? + 62 Ie) exp(6), > 1,61, 62,0 > 0,2 > 0, 


where 6}, 62, 8, C» are independent of n. 


4 Convergence 


In this section, we study the convergence of the proposed Crank—Nicolson 
scheme for (1) with initial and boundary conditions (2)—(4). 


Let u?; (1 <i<I-1,1<j<J—1,n=1,2.,...) be the exact solution of 


(17) and (18), and let UP, (1<i<I-1,1<j<J-1,n=1,2,...) be the 
exact solution of (15) and (16). Define €?; = Uj; —uz,; (l<i<I-1,1< 
j<J—1,n=1,2,...). By considering €" instead of €?, we obtain 
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Agi? 1 ¢! Ai? 1 1 
( ) x = gE ( ) x = él él 
T(3— 6) 22-8 At [(3— 6) 22-8 2 (34) 
1 
= + O(Ax)? + (O(Ay)? + O(At)?-*), 
and 
(anime ff ate at pene 
T3- BY ee nk) —hy (ga=a — 8G 
‘i Em aenet 
92-8 Ag 
(At)i-8 fe 1 a 1.) (35) 
+TG—p) {- a (bn—k—1 — bn—n) &* — (33-8 — bert + yenaé 
aS en-l 4 En 
) 
- AEr—1 + AEr 


. + (O(Az)? + O(Ay)? + O(At)?-8), n> 2. 


Now, we are ready to present our next theorem. 


Theorem 6. If €* € Hé(Q), then the solutions of the finite difference ap- 
proaches (17) and (18) are unconditionally convergent. 


Proof. Let a = _ Ay" If we multiply (34) by (or) then we obtain 


as ees Pe 1 eas 
ean eee 
a (O(Az)? + O(Ay)? + O(At)?-*) S50. 


(36) 


Applying Theorem 2 (Green’s theorem) to (Aé',€+) in the left side of (36) 
and applying Theorem 1 (the Cauchy—Schwarz inequality) and Lemma 2 to 
the right side of (36), we find that 


lee I 


1 
arrap lt I + allel + + 5llverll’ 


_ lox? ay +o(aeyFl jet 
ES 9 : 


Therefore, 
IVE]? < ]O(Ax)? + O(Ay)? + O(At)?-? |? (37) 


If we multiply (35) by (€" — €"~1), then we obtain 
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a 


PFA (é" _ En-1 en _ oy ai a ee _ co 


4 ; (é” 4 eek ge ame ; (AEé” Ag”! €” _ or) 
ee a ae 


(38) 


‘Gee - th) — ae) 
n—2 

+a Ne (by—K-1 = bn—k) ae 3s oe 
k=1 


+ ogra = b1) (ence Ome et} 
+ < (O(Az)? + O(Ay)? + O(At)?~*) ,€" — E71). 


Again, using Theorem 1 (the Cauchy—Schwarz inequality) and Lemma 2, we 
can write 


< (O(Ax)?+O(Ay)? + O(At)?-9), €? — €™1 > 


cael 
lg , 


n||2 
< |O(Az)? + O(Ay)? + O(At)?-9 ||? "I ; 


Simplifying relation (38) (similar to Theorem 5, in which the simplification 
of (23) resulted in (29)) and using the recent relation, we obtain 
1 a 


2” 2-8 


en? +5 Ive 


2 n—2 
< |]O(Ac)? + O(Ay)? + O(A#)-? IP? + SS neat — One) [le 
k=1 


n—2 
ay? 2 aD ay? fs ened 
+ Tae De naka — bre) (TEA + MEE) + ae lle? 
k=1 
ary? ay? a soda iD on decth 
(> + pet aoe TONE TW +a ilver "|, 222. 


By Theorem 3 (the Poincaré-Friedrich inequality), there exists a constant 
Cp > 0 such that 
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1 a 
5 IVEr||? < ]O(Ax)? + O(Ay)? + O(At)?-F |? 
2 n—2 


+ —=~Cy (bn—k—1 — bn—x) || VEF|I? 
k=1 


n—2 


ay? ky2 b-1112 
a (Aspe? 2 nko ~ bre) (IVE I? + || ver)?) a 


ep || Ve" 


ay? 


Bip 
2D 2 

ay ay a 

(> + Tape t 2x28 


lve" 


t Lep || Ver} |]? 


2 
eo EBED, 


1 
2 
As we know, €° = 0. Without loss of generality, relations (37) and (39) can 
be written as 
||VEt||” < ]O(Aw)? + O(Ay)? + O(At)?-F |’, 
n—-1 
II VE"||? < 2|]O(Ae)? + O(Ay)? + O(At)?-F ||” + S> Cy || VER’, 
k=1 
n>2, Ch>0 for k=1,...,n—1. 


(40) 


Thus, by using Theorem 4 (the discrete Gronwall theorem), the set of equa- 
tions (40) yields 


n-1 
IVEr||? <2 ||O(Az)? + O(Ay)? + o(At)?-8 |) exp & c,) , nol, 


k=1 


and according to Theorem 3 (the Poincaré—-Friedrich inequality), there exists 
a constant € > 0 such that 


n-1 
ll"? < 22, ]O(Aa)? + O(Ay)? + O(At)?-F ||? exp & c:) « eat, 


k=1 


(41) 
By using Lemma 3, it is easy to show that 
= Say? a 
2 


Set ¢ = (207 + ae aoe 2) Cp+1. Then, relations (41) and (42) allow 
us to write 
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IE" Il < \/ 2p exp(¢) ||O(Ax)? + O(Ay)? + O(At)?-* |], n>1,¢>0,%)>0. (43) 


5 Numerical experiments 


In this section, we present some numerical tests that confirm the validity of 
the proposed numerical method. To measure the accuracy of the proposed 
method, we use the maximum absolute error given by 


Loo = max U;3(T) a U;i,5(T) , 


1<i<E1<j<J 


where U;,;(T) and U;,;(T) denote the numerical solution and the exact solu- 
tion of (1) with initial and boundary conditions (2)—(4) at (x, y;) and time 
T, respectively. 


Example 1. Consider a two-dimensional test problem of the form (1), 


with 2 = [0,1] x [0,1], f(#,t) = (BES + 2h. + 2ttn?) sin(mz + my) + 


t* sin(r2 + my), and suppose that the initial and boundary conditions are 
assumed using the exact solution u(z, y,t) = t*sin(x + my); see [13]. Now, 
we provide some tests. 


Test 1 Kumar, Bhardwaj, and Dubey [13] considered this example using 
a local meshless method with 2025 points on 2. They reported the maximum 
absolute errors and CPU time with 8 = 1.7,1.9 and different values for At 
at the time 7’ = 1.0. Using the proposed method, we repeated this test. We 
considered this example by assuming J = J = 45 (2025 points on Q). To 
solve the linear system of equations, we used the GMRES-m method with 
m = 20. 

Table 1 presents the maximum absolute errors and CPU time obtained by 
Kumar, Bhardwaj, and Dubey [13] and the results of the proposed method 
with 8 = 1.7, different values for At, and 2025 points on [0,1] x [0,1] at 
T = 1.0. Table 2 presents the maximum absolute errors and CPU time 
obtained in [13] and the results of the proposed method with 6 = 1.9, different 
values for At, and 2025 points on [0,1] x [0,1] at T = 1.0. 

As Tables 1 and 2 show, the maximum absolute errors and the CPU time 
of [13] and those of the proposed method are close, but the CPU time of the 
proposed method is smaller than that of [13]. 

The following tests show that the proposed method provides acceptable 
accuracy with a smaller number of points on (Q. 

Test 2 We considered this example by the proposed method with Ax = 
Ay = 0.1, 6 = 1.5,1.9, and different values for At. According to Table 3, 
with different values for At, the maximum absolute errors were small enough 
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at T = 1.0. Also, decreasing the size of the time step increased the CPU time 
very slowly and improved the accuracy. The value At = a was selected for 
the next test. 

Test 3 We considered this example by the proposed method with At = a 
B = 1.5,1.9, and different values for Ax, Ay. According to Table 4, the 
accuracy was acceptable. Also, the CPU time was reasonable with Ar = 
Ay = un 0° Moreover, by decreasing Ax and Ay to a 30° the CPU time 
increased rapidly, and the accuracy did not improve significantly. According 
to relation (43), the convergence rate of our method depends on Az, Ay, 
and At. In this case, the space steps decrease, but the time step is constant. 
Therefore the accuracy does not improve. 

As shown in Tests 2 and 3, a very small size the of space step is not 
recommended, but small size of a time step is recommended. According to 
Diethelm, Garrappa, and Stynes [7], a high-order space discretization for a 
time-fractional partial differential equation is not advisable. They believe 
that to reach a high convergence, we must choose very small size of the time 
step in comparison with the size of the space step. Our experiments confirmed 
this idea. 


Table 1: Comparison of the maximum absolute errors and CPU time with 
(@ = 1.7, different values for At, and 2025 points on [0,1] x [0,1] at T = 1.0 


At] Lo [13] om GPU (s) [13] CPUs) 
as 1.2917e —02 1.2575e — 02 1.751 1.414 
a 5.4532e —03 8.7100e — 03 2.210 1.996 
aa 2.335le —03 5.0648e — 03 3.062 2.746 


Table 2: Comparison of the maximum absolute errors and CPU time with 
(G =1.9, different values for At, and 2025 points on [0,1] x [0,1] at T = 1.0 


At] Le [13] i GPU (s) [13] CPUs) 
a 2.7619e —02 1.6456e — 02 1.751 1.298 
wi 1.3079e —02 1.0294e — 02 2.210 1.613 
a 6.19538e — 03 5.7166e — 03 3.062 2.119 


6 Conclusion 


The Crank—Nicolson difference scheme can be used easily for space-fractional 
equations, but some manipulations are needed for time-fractional equations. 
In this paper, the Crank—Nicolson method was extended for the discretiza- 
tion of a TFTE. The solvability, stability, and convergence of this proposed 
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Table 3: Maximum absolute errors and CPU time for different values of At 
and @, with Ax = Ay = 0.1 at T = 1.0 


BH=15 B=19 
At Te CPU(s) Tee CPU (s) 
4 | 1.3159e—02 0.1069 | 1.9190e—03 0.1041 
+, | 1.0567e—-02 0.1249 | 1.2971le—03 0.1269 
d+ | 7.500le—03 0.1673 | 8.3675e-03 0.1713 
+, | 5.5027e—03 0.2944 | 5.7025e-03 0.2908 
790 4.3813e —03 0.6513 | 1.9956e—03 0.6423 
aa | 3-7898e 03 1.7490 | 3.3609e—03 1.7547 


Table 4: Maximum absolute errors and CPU time for different values of Az, 
Ay, and 8, with At = ry at T= 1.0 


B=15 B=19 
Aa = Ay Loo CPU(s) Lig CPU (s) 
- 5.5027e — 03 (0.2944 | 5.7025e —03 (0.2908 
aa 3.1200e—03 0.4722 | 3.6015e—03 0.4273 
a 2.5315e— 03 3.6787 | 3.0738e—03 = 2.2135 
a 2.3744e —03 74.7221 | 2.9432e—03 35.9883 


method were proved. The numerical results were accurate enough. Accord- 
ing to the numerical tests, to reach a high convergence, a very small size 
of the space step is not recommended, but a small size of the time step is 


recommended. 
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